Anisotropic anomalous diffusion modulated by log-periodic oscillations 
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We introduce finite ramified self-affine substrates in two dimensions with a set of appropriate 
hopping rates between nearest-neighbor sites, where the diffusion of a single random walk presents 
an anomalous anisotropic behavior modulated by log-periodic oscillations. The anisotropy is revealed 
by two different random walk exponents, i/x and Vy, in the x and y direction, respectively. The values 
of these exponents, as well as the period of the oscillation, are analytically obtained and confirmed 
by Monte Carlo simulations. 
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I. INTRODUCTION 

The underlying mechanisms of anomalous diffusion on 
fractal structures has attracted the attention of scien- 
tists for many years (see, for example Ref. [T] and ref- 
erences therein). In this regard, it has been recently 
found that, on some kind of self-similar substrates, in ad- 
dition to the well-known subdiffusive behavior, the mean- 
square displacement of a random walk (RW) is mod- 
ulated by logarithmic periodic oscillations [2-4J. The 
same kind of modulation was also observed in biased 
diffusion on random systems ^, earthquake dynam- 
ics [5] , escape probabilities in chaotic maps [7] , processes 
on random quenched and fractal media [8], diffusion- 
limited aggregates [S] , growth models [TU] , and stock mar- 
kets [TT]. There is general agreement that this ubiqui- 
tous phenomenon appears because of an inherent self- 
similarity ^12], responsible for a discrete scale invari- 
ance |T3j. Nevertheless, this self-similarity has to be iden- 
tified for every system. 

The origin of log-periodic modulation can be easily de- 
termined for a minimal model of RW introduced in [3]. 
This model, which depends on two parameters, L £ N 
and < (5 e M, consists of a one-dimensional lattice 
and a single particle moving by jumps between nearest- 
neighbor (NN) sites. The hopping rates are defined in a 
way that a region of size L" (with n = 0, 1, 2...) is char- 
acterized by a diffusion coefficient D^"), and the ratio 
between any two consecutive coefficients is a constant, 
i. e., £)("+i)/D(") = S for all n e N. As a result, the RW 
mean-square displacement is modulated by log-periodic 
oscillations and, both the RW exponent and the period 
of the oscillations can be obtained using rather simple 
arguments and calculations (for more details, see [3]). 

This method can also be applied to the study of RW 
on a self-similar substrate in two dimensions. It has been 
shown [_4 that, in this case, each region of size L" x L" (L 
is the basic length of the substrate, and n = 0, 1, 2, ...) is 
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characterized by a diffusion coefficient £)("). Here again 
a subdiffusive behavior modulated by log-periodic oscil- 
lations arises, because the ratio Z3("+i) /£)(") takes a con- 
stant value. It is the symmetry between x and y direc- 
tions which allows the heuristic arguments used in the 
one-dimensional case to be easily generalized to calcu- 
late de values of the RW exponent and the period of the 
oscillations. The important point is that, for a particle 
in a central square of size L" x L", the typical time to 
leave this square along the x direction is the same as that 
along the y direction. 



In this paper we investigate single particle diffusion 
on self-affine structures. In general, the lack of symme- 
try between the two main directions (x and y) makes the 
analytical treatment difficult. However, the problem sim- 
plifies considerably for a special kind of substrate, that 
in which the space explored by a RW grows with the 
same anisotropy as the substrate itself does. We study 
this case first. The same kind of arguments employed 
to analyze diffusion on self-similar substrates allows us 
to show that, in this case, the mean-square displacement 
as a function of time is a power-law modulated by log- 
periodic oscillations but, in contrast with its self-similar 
analog, the specific properties of this function are now 
direction-dependent. Indeed, although the period of the 
modulation is isotropic, two different RW exponents ex- 
ist, one for the displacement in the x direction, another 
for the displacement in the y direction. We compute an- 
alytically the RW exponents and the period of the mod- 
ulating oscillation, and confirm these results by Monte 
Carlo simulations. 



For the sake of completeness, we then study numeri- 
cally the RW behavior on a more general self-affine sub- 
strate. The outcomes of these simulations suggest that, 
also here, the mean-square displacements along the x and 
y directions, as a function of time, follow log-periodic 
modulated power-laws, which are independent of each 
other. 



II. ANALYTICAL APPROACH 

We study the behavior of a RW on two self-afRne sub- 
strates, referred in what follows as model I and model 
II. Each substrate is built in stages, and the result of 
every stage is called a generation: a periodic array of 
basic or unit cells which consists of sites connected by 
bonds. We denote by L^ and Ly the linear size of the 
unit cell of the first generation in the x and y directions, 
respectively. On these substrates the motion of a single 
particle occurs stochastically. At every time step, the 
particle jumps with a non-zero probability only between 
NN sites which are connected by a bond. The details of 
each models are given below. 



A. Model I 

The building process is illustrated in Fig. [l] which 
shows the unit cell for the zeroth, first, and second gen- 
eration. It is easy to see that, for this model, L^, = 5 and 
Ly = 3, where the length unit is the distance between 
NN sites. It is also apparent from this figure that the 
second-generation unit cell has linear sizes L^ and L^, in 
X and y direction, respectively, and is built from the first- 
generation one in a self-affine way. In general, the linear 
sizes of the nth-generation unit cell are L" and L" and 
the corresponding two-dimensional periodic substrate is 
obtained by connecting these cells (the first-generation 
substrate is sketched at the top of Fig. [2]) 



The hopping rate between any NN connected sites in 
the X direction is always fc^*^'. On the other hand, the 
hopping rate in the y direction depends on the site and 
on the generation. Their values are determined by ask- 
ing that the mean time to leave a nth-generation unit 
cell along the x and y directions coincide. We call i*-"^ 
this escape time. Because of this constraint, there will be 
n + 1 different hopping rates (fc^'^, i = 0, ..., n) related 
to the nth generation. As an example, in Fig. [l]we show 
an schematics of the the zeroth, first and second gener- 
ation, with one, two and three kinds of hopping rates, 
respectively. In this sketch, a thin bond represents k^^\ 
while the other hopping rates are represented by thicker 
bonds. We can observe that k'^^\ appears at the top of 
the first generation unit cell, and k^^i appears at the top 
of the second generation one. 

We proceed now to analyze the behavior of the diffus- 
ing particle on a nth-generation substrate. It is useful 
to remember that, on any periodic substrate, normal dif- 
fusion should be observed if time is long enough for the 
RW to be influenced by the structure periodicity. As 
we work with an asymmetric substrate (i.e L" ^ L" 
for the nth-generation) we have to consider x direction 
and y direction separately. For the nth-generation sub- 
strate, a diffusion coefficient Dx {Dy ) in x (y) di- 
rection can be defined through the time dependence of 
the mean-square displacement A^a;(t) = {[x{t) — x{Q)Y) 
(A'^yit) = {[y{t) - 2/(0)]2)), i.e., via the relations 



A^x{t) - 2i?(")i, 



(1) 



and 




FIG. 1. The unit cells of model I. The zeroth, first, and second 
generation are drawn in (a), (b), and (c), respectively. The 
basic length-scales are Lx = 5 and Ly — 3. A thin bond 
(thick bond) represents a hopping rate fc'"' {k^'' , i > 1). 
More details in the text. 



The full self-affine substrate, we are interested in, is 
the result of an infinite number of iterations. Note that 
this substrate is finitely ramified, and that a region of 
size L" X L" can be separated from the rest by cutting 
four bonds. 



A'y{t) = 24")i, 



(2) 



valid for a time t longer than t^"' . 

The diffusion problem is trivial on the zeroth- 
generation substrate. This is a simple square lattice, and 



J, y 



(3) 



The first-generation substrate (top of fig. [2]) presents a 
more difficult task. However, regarding x-direction diffu- 
sion, the whole substrate and the string of cells displayed 
at the bottom of the same figure, with periodic boundary 
conditions in the y direction, lead to equivalent problems. 
We exploit this equivalence and calculate the diffusion 
coefficient of that one-dimensional array, following the 
steady-state method [15j. We get 



J^ir 



and thus, 



^i"V^i" 



+1) 



(-)"fc(°\ forn = 0,1,2,..., (4) 



7/5, forn = 0, 1,2,... . (5) 




Periodic Boundary Conditions 



i 



i 



FIG. 2. First generation of model I. Top: the substrate built 
with the basic cell shown in Figlllb). Bottom: the infinite 
one-dimensional string of cells used to compute the diffusion 
coefficient Dy . The arrows indicate periodic boundary con- 
dition in the y direction. For example, if a RW at site a 
(b) jumps upward (downward), with a hopping rate fc'^', it 
arrives at site b (a). 



To find the diffusion coefRcients in the case of the y 
direction, we divide Eq. (pi) by Eq. (fTl), imposing the 
same escape time constraint, i.e., A^a;(t'"-') — L^" and 



A^yit^'^')) = L2". This leads to 



D 



in) 



T2n 
_ ^ 

where the diffusion coefficients 

4") = (^rfc(0), for n = 0,1, 2, 



(6) 



(7) 



can be obtained from (using Eq. (B and the values of 
Lx and Ly). Hence, the ratio between consecutive coef- 
ficients is also a constant in the y direction: 



DM/Dij^+i)^Sy=35/9, for n = 0,1, 2, 



(8) 



At this stage, the model is completely defined, and the 
hopping rates are obtained recursively from ([5]), by using 
the above mentioned trick of converting the diffusion two- 
dimensional problem in a one-dimensional problem: 
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FIG. 3. (Color online) Schematic of the mean-square dis- 
placement in the x direction, as a function of the time, shown 
by the red (tick) curve. The length of the segment be is 
logio(2Di"') -logio(2Di"+')) = logio(5,), because of Eq.^. 
From the slopes (= 1) of the full straight lines (represent- 
ing the normal diffusion behaviors, A^x = 2Di" i), one gets 
that the segments ad and cd have the same length or, equiv- 
alently, that logj^Q(T) = logjQ(L^) -I- \q%^q[&^"). The dashed 
straight line represents the global power law I^x ~ t^"^ with 
2v^ = logigL^/logio-T. Thus, v^ = (2 + logi(, (J^/logjo L:i,)~\ 
The mean-square displacement in the y direction exhibits an 
analogous behavior. 



Let us now consider a RW on the full self-affine struc- 
ture. For a time t in the interval [i^"), t^""*"^^], the follow- 
ing relations hold 



L-<./W^)<Vl^\ 



L-<./Wy(t)<Ll 



(10) 



(11) 



(9) 



and it will be impossible for the RW to distinguish the full 
self-affine structure from the nth-generation one. Thus, 
Eqs. (IT]) and (pi) account for the RW behavior in that 
time window, and the mean-square displacement should 
behave qualitatively as sketched in Fig. [3] This behavior 
is reminiscent of single particle diffusion on a self-similar 
substrate, whose mean-square displacement as a function 
of time obeys a log-periodic modulated power-law [3]. 
Because of the lack of symmetry between the x and the 
y directions, to describe diffusive behavior in the case of 
a self-affine substrate, we need not one but two functions, 
which we expect to be 



and 



A'x{t) = C,t"'-Mt) 



A'y{t) = Cyt'-'y fy{t) 



(12) 



(13) 



where C'x and Cy are constants, Vx and Vy are the RW 
exponents, and fx{t) {fy{t)) is a log-periodic function 
with period Tx (Ty). 

The values of these quantities can be computed from 
the parameters of the model, after simple geometrical 
analysis of Fig. [s] (see figure caption and Refs. [3j 1^ for 
further details). The results are 



1 

log 10 S^ 

logio ix 



1 

l°Sio ^y 
logio Ly 



(14) 



(15) 



Model II 



For this model, the unit cells for the zeroth, first and 
second generation are shown in Fig. [4] The full self-afHne 
substrate is here also obtained when the generation order 
goes to infinity. The linear sizes of the rath-generation 
unit cell are L" and L" with Lx = 3 and Ly — 2. 



(a) 



(b) 



r 



, (1) 



r 



(c) 

.(2) 



T^^ 



and 



Tx — OxL^ 



Ty — ^y^y 



(16) 



(17) 



Note that, even when Vx 7^ ^y, the period of the mod- 
ulations coincide, because of the constraint ^, i.e., 



= S.Li = 



D 



(n) 



(n+1) 



Ll = 
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D 
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where we have also used ([5|) and (pi). We call r this 
period. From the equations above, the values of the pe- 
riod and the exponents are r = 35, Vx = 0.4527 and 
i^y = 0.3090. 

Let us note that the average time to escape from a unit 
cell of the nth-generation is i*-"' = r", which means that 
the relations (10) and (11) hold for. 



<i< 



^n+l 



(19) 



Then, when the RW leaves the initial region, of size 
L" X Ll^,, to enter the next one, of size L""*"^ x £"+^, the 
length-width ratio {L^/Ul) is increased by an anisotropic 
factor a = Lx/ Ly, while the average time increases from t 
to ri. On the other hand, according to ( 12 1 and (14), the 



corresponding mean square displacements are related by 
A2a;(rt) = LlA^x{t) and A^y{Tt) = LlA^y{t). There- 
fore, in this transition, the ratio y^A'^x/A-^y is also in- 
creased by a factor a; i.e., the space explored by the RW 
grows with the same anisotropy as the substrate where 
the diffusion takes place. 



FIG. 4. The unit cells of model II. The zeroth, first, and 
second generations are shown in (a), (b), and (c), respectively. 
I/a; = 3 and Ly — 2. 

The diffusion of a single particle is analyzed as on 
model L That is, we reformulate the two-dimensional RW 
problem on a one-dimensional array and compute the dif- 
fusion coefhcients following the steady-state method [IS] . 

For the nth generation, we obtain 



^(n) = (_)«fc(o)^ for n = 0,1,2,..., (20) 

and thus, 



£'(,"V^i"+^' =Sx^ 4/3, for n = 0, 1, 2, ... 



(21) 



In average, the time to leave a nth-generation unit cell 
along the x direction becomes the same as that along the 
y direction if 



Din) 






for n = 0,1,2,... 



which implies 



i?("V4"^'^ = '5^ = 3, for n = 0,1, 2,. 



(22) 



(23) 



The fc'*''s, coming from (22 1, are again computed from 
([9]) (with Lx = i and Ly ==2). Furthermore, in spite of 



the differences between model I and model II, the quali- 
tative behavior sketched in Fig. [3] we expect to be valid 
for both models. Therefore, the RW exponents v^ and 
Vy, and the period r are given by ([l4|, (fTsj) and ([Tsl); 
with the values ; 



0.4421, 



0.2789 aSd 



12. 



III. 



NUMERICAL RESULTS 



To test the predictions outlined above, we perform 
standard RW Monte Carlo simulations, on a nth- 
generation unit cell for each model. In model I (II) every 
RW starts at the center of symmetry of the cell (at the 
top- left most site) . The value of n is always chosen large 
enough to prevent the RWs from reaching the cell borders 
(the bottom and right cell borders) during the simulation. 
Working on this cell is thus equivalent to working with 
the full self-afBne structure. In all simulations the hop- 
ping rate k^^' is set to 1/4, and the other k^^'^s (i > 1) 
are obtained from ^. After every Monte Carlo step, the 
time is increased by At = 1 

With the numerical results of model I, in Fig. [5]-(a) 
we have plotted the mean-square displacement along the 
main directions. We see in these plots that both A'^x{t) 
and A^y{t) are well described by modulated power laws. 
The upper and lower straight lines have slopes 2z^2; and 
2i'y, respectively. They are drawn to guide the eyes, us- 
ing the analytical values of the RW exponents. The log- 
periodicity of the modulations can be better observed 
in Fig. [5]-(b), where we have plotted \ogig{A'^x/Axt'^'^''') 
and \ogiQ{A'^ y/Ayt'^'^y) against log;^o(^)j using the same 
data as in the part (a). Ar^ (Aj^) is a constant cho- 
sen to have the oscillations in the x{y) direction cen- 
tered around 0.00 (0.05). The continuous lines are of 

the form S sin(27r log]^Q(i)/ log]^o(''') + Q^)' i-^-, the first- 
harmonic approximation of a periodic function with pe- 
riod log;^Q(T), where B and a are fitted parameters and 
T = 35 (the above given analytical period). It is clear 



from this figure that the theoretical predictions (Eq.(14|, 
( 15 1 and p^) are consistent with the numerical findings. 



The corresponding numerical results for model II are 
shown in Fig. [6] Note that, also for this model, at 
long times, the mean-square displacement as a func- 
tion of time is well described by modulated power 
laws. To better appreciate the log-periodicity of 
the modulation, we have plotted \ogiQ{A'^x/Axt'^'^^) 
vs. logio(0 aiid logiQ^A'^ y/Ayt'^'^y) vs logio(0 in the in- 
set of this figure. The fitting curves are of the form 

Dsin(27rlog]^o(^)/logio(''')+'^)' with the analytical value 
T = 12. The agreement between analytical and numerical 
results is also good. 

We consider now a substrate (model III) which con- 
sists of the full self-affine structure of model I but with 
the same hopping rate k^^' between any pair of con- 
nected NN sites. For this model, the average time to 
leave a n-generation unit cell along the x direction is dif- 
ferent from that along the y direction. It may occur that 
LS < yA%(i) < LS+\ L^ < ./^W) < L^^\ for a 




0.1 



X 



,.0.0 



-0.1 



- □ X 




J^^ rrtBffi 




^yx^ 


\/\^ 


O ^J^ 


\/\/ ' 


o 




CP 




o 


(b) 


- 


1 1 1 



10" 



10 



10' 



10° 



FIG. 5. (Color online) (a): The mean-square displacements 
A^x (green squares), and A^y (purple circles) as functions of 
time for Model 1. The upper strai ght line has a slope 21/^ , with 
Vx = 0.4527 obtained from Eq. (141. The lower straight line 
has a slope 2i'y, with Vy — 0.3090 obtained from Eq. ( 15 1, (b): 



X = log.oiA^x/Axt^-'^) vs. logiot (y = log,o(A2y7l',f2-„) 
vs logjQ(i)) for the same data. Ax and Ay are a properly 
chosen constants. The curves represent the first-harmonic 
approximations Ba; sin[27rlogj^Q(f)/logjQ(r)] +a (blue-upper) 
and By sin[27rlog]^)/logio(r)] +13 (red-lower). The period 
T is given by Eq.(|18[|. Bx, By, a and /3 are fitted constants. 



given time t and m ^ n, which, in other words means 
that, near i, the RW behaves as in the rt-generation 
substrate, regarding the x direction, but as in the m- 
generation substrate, regarding the y direction. Thus, 
we cannot expect the heuristic arguments in the previous 
section continue to be valid and we have then to study 
the problem numerically. 

The logarithm of the scaled mean-squared displace- 
ments (in the x and y directions), i.e., A'^x/{Axt'^'''') and 
A^y/{Ayt^''y) are plotted in Fig. [t] as a function of the 
logarithm of time. The RW exponents v^ = 0.4373 and 
Vy = 0.3859 in this figure are fitted values . Let us note 
that Vj. is different from Vy , and that the data of Fig. [71 
strongly suggest that the modulation have the same pe- 



y direction Vy, and that the global subdiffusive behavior 




FIG. 6. (Color online) The mean-square displacement A^x 
(black squares) [A^y (red circles)] versus time for model II. 
The top straight line has a slope 2ux ~ 0.8842, and the lower 
straight line has a slo pe 2i/ y = 0.5579. Both exponents are 
obtained from Eq.([l4| 
logio ^^x/A^t^"- 
(red circles) against t, for the same data. The curves were ob 
tained as in Fig.|5| with the period t calculated from Eq.( 18 1 



and (151 



The inset are plots of X = 
(black squares) and Y = log^g A^y/A^f^'''" 



riod r in both directions. As expected, the numerical 
values of these parameters are not in agreement with 
Eqs.([l4]), ([15]), and ( |l8t. W e would like to remark that if 
we used Eqs.([l4| and |l5[) (with S^ = 7/5 and 6y = 7/3 
resulting from the new hopping rates) we would get the 
RW exponents v'^ = 0.4527 and ly'y = 0.3609, which, in 

turn would lead to the periods r^ — LJ ^ = 35 and 
Ty = Ly " = 21; different of each other (see caption of 
Fig. ^for the equation r = L^^"). Note that the numeri- 
cal value of i^x ~ 0.4373 {vy — 0.3859) is smaller than i/'^ 
(larger than i/'), and the numerical value period t is in 



the range \tLt^] (r = 26). For model III, the Eqs. (10 1, 
(11) and (19) do not hold because, in average, the KW 



reaches the top or bottom border of the n-generation unit 
cell before reaching the right or left border of the same 
cell. In the case of model I, this is avoided by properly 
modifying some hopping rates in every generation. The 
diffusion spread in the y direction is thus slowed down 
(A:^"^ < fc("~^\ see Eq.(|9l), and the horizontal and verti- 
cal cell borders are, in average, simultaneously reached. 



IV. 



CONCLUSIONS AND DISCUSSION 



We have studied the problem of single particle diffusion 
on a finitely ramified self-affine structure in two dimen- 
sions. For a special kind of models, for which the ratio 
between the x and y mean-square displacements matches 
the structure anisotropy, we argue that the RW expo- 
nent in the x direction v^ is different from that in the 
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FIG. 7. (Color online) Scaled mean-square displacements for 
Model III. Plot of X = logio A^x/Caii^"" vs. log^o < (green 
squares) and Y — log^g A^y/Cyt^"^ vs. logj^Q t (orange cir- 
cles), using numerical data {Cx and Cy are properly cho- 
sen constants, see the text). The full lines represent the 
first-harmonic approximations A3; sin(27rlogioi)/logio(''") -I- 
a) (upper) and Ay sin(27rlog^o i)/logio (''")+ /3) (lower) of X 
and Y, respectively. Here, t = 26, Ax, Ay, a and /3 are 
fitted constants. The uppers dashed line represents the first- 
harmonic approximation of X, with period r^ = 35. The 
lower dashed line represents the first-harmonic approximation 
of Y , with period Ty = 21 



is modulated by log-periodic oscillations with a period t 
which does not depend on the direction. The arguments 
employed in this work allow the main properties of the 
particle mean-square displacement to be obtained as a 
function of model parameters. Because our arguments 
are somehow heuristic, MC simulations using two mod- 
els, I and II, were also carried out. The numerical results 
confirm our theoretical predictions. 

For the rest of the self-similar systems, our conclu- 
sions are more limited, due to the lack of suitable an- 
alytical methods and that the RW explores the space 
with an anisotropy different from that of the substrate. 
The results of the MC simulations performed using one 
of these models (III), show (within the accuracy of the 
simulation) that, also in this case v^ ^ Vy and the RW 
mean-square displacement is modulated by log-periodic 
oscillations with an isotropic period. However, we can- 
not guarantee that this behavior will hold in the limit of 
an arbitrary long time; that is why we have introduced 
models I and II. Let us finally note that the extension 
of our analytical results to other values of L^ and Ly is 
straightforward. 
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